Mon. Not. R. Astron. Soc. 000, [|-| (2000) Printed 1 February 2008 (MN WI^ style file vl.4) 



Flaring and warping of the Milky Way disk: not only in 
the gas.... 
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ABSTRACT 

This paper presents an investigation of the outer disk structure by using data from 
a recent release of the 2 micron sky survey (2MASS). This 2MASS data show unam- 
biguously that the stellar disk thickens with increasing distance from the Sun. In one 
of the field (longitude 1=240) there is also strong evidence of an assymetry associated 
with the Galactic warp. This flaring and warping of the stellar disk is very similar 
to the features observed in the HI disk. The thickening of the stellar disk explains 
the drop in density observed near the Galactic plane: stars located at lower Galactic 
latitudes are re-distributed to higher latitudes. It is no longer necessary to introduce 
a disk cut-off to explain the drop in density. It is also not clear how this flaring disk 
is distinct from the thick disk in the outer disk region. At least, for lines of sight in 
the direction of the outer disk, the thickening of the disk is sufficient to account for 
the excess in star counts attributed by some models to the thick disk. 

Key virords: Galactic structure 



1 INTRODUCTION. 

The structure of the stellar disk of our Galaxy is stiU not 
very well known. The stellar population of the disk is usu- 
ally decomposed in a thin disk and thick disk component. 
The thin disk and thick disk populations are separated on 
the basis that they may have different metallicity and ve- 
locity profiles. The typical ratio of the thick disk to thin 
disk components in the solar neighborhood is about a few 
percent (Gilmore & Reid 1983). . Both disk are supposed to 
follow a double exponential profile, associated with a radial 
scale length and a vertical scale height. It is unclear at the 
moment if the scale length of these 2 disks are the same. 
There is still some uncertainty even for the scale length of 
the thin disk itself, although the more recent measurement 
seems to indicate that it is between 2 and 3 Kpc. For the the 
thick disk the uncertainty is much larger, typically the scale 
length is in the range 1.5 to 4.5 Kpc. The scale heights of 
the thin disk is about 0.3 Kpc , and is around 1 Kpc for the 
thick disk. The structure of the gas component (whether HI 
or CO) is better understood, and suggests some questions 
concerning the stellar disk. Since a warp is well identified in 
the gas ( W. B. Burton and P. te Linkel Hekkert 1986), one 
may wonder if it does exists also in the stellar components, 
with the same amplitude and phase. Another problem is the 
extension of the gas component which is visible up to a dis- 
tance of 25 Kpc from the Galactic center. Some observations 
suggests that the stellar disk ends long before, at about 14 
Kpc from the center. This fact may indicate that stellar for- 
mation in the disk did not occur beyond 14 Kpc, or that 



Table 1. The 3 fields 



Field 


field center 


size of field 


Number of stars 


L66 


(l,b)={66,0) 


(Al,Ab)=(2,100) 


1.4 lO'' 


L180 


{l,b)={180,0) 


(Al,Ab)=(10,100) 


3.8 10" 


L240 


{l,b)={240,0) 


(Al,Ab)=(10,100) 


4.7 lO** 



the geometry of the stellar disk is not well understood. One 
additional source of complication is the possibility that the 
disk in addition to its warp mode presents an elliptical mode 
of deformation. At the moment this issue is very speculative, 
however this effect may exist and should be investigated. 



2 THE DATA 
2.1 Presentation 

The data set consist of 3 strips perpendicular to 
the Galactic plane around the longitudes, 1=66, 1=180, 
1=240. AH the data have been extracted form the 
2MASS point source catalogue acces sible on the Web 
(tittp:/ /www. ipac.caItech.edu/2MASS/). The detailed char- 
acteristics of the fields are given in Table 1. 
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Figure 1. Mean extinction in the V band in the 3 fields. These 
plot were constructed from the map of Schlegel et al. (1998). 

2.2 Canceling the effect of extinction. 

By combining 2 photometric bands it is easy to derive a 
magnitude that is independent of reddening. For instance 
using J and K, we can derive Me- 

Me ^K- , ^"^^ (J - K) 
A.J - Ak 

Following Rieke and Lebofsky (1985) we take: 

-P- = 0.112 and — ^ = 0.282 

Av Av 

Leading finally to: 

Me = K - 0.659 (J - K) (1) 

This corrected magnitude will be very insensitive to redden- 
ing, furthermore, it is important to notice that the reddening 
in the 3 fields is low (see Fig. 1). The maximum reddening 
is about Av — 1.5, which is only 0.17 in the K band. Conse- 
quently, after adding the correcting term: —0.659 (J — K), 
the residual will not exceed a few percent. 

2.3 Star counts as a function of latitude 

Before making the star counts, it is essential to correct for 
empty areas in the catalogue (due most of the time to the 
obliteration of fainter stars by a very bright star). To locate 
these bad areas a 2 dimensional histogram in the 1 and b co- 
ordinates was computed. The size of the bin was chosen in 
order to have a mean of 50 stars in the bin. It was required 
that the bin contains a minimum of 5 stars to be considered 
valid. This way the histogram of the valid areas was con- 
structed as a function of latitude. The step in latitude is 0.5 
degree. The counts were finally normalized to a surface of 
1 square degree. To avoid any problem related to the com- 
pleteness of the sample, only stars brighter than Me = 12.5 
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Figure 2. Star counts for stars brighter than Me = 12.5 as a 
function of latitude for the 3 fields. The counts are normalized in 
number of stars per square degree. 
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Figure 3. Difference (percentage) between counts at positive 
and negative latitude as a function of latitude. 

have been included. The resulting histograms are presented 
in Fig. 2. 
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2.4 Assymetry between positive and negative 
longitude. 

It is interesting to investigate the symmetry of the distri- 
bution of the counts as a function of latitude. A significant 
assymetry would reveal the presence of the Galactic warp. 
In Fig. 3 the difference in star counts between positive and 
negative longitude is plotted. It appears immediately that 
the field L240 presents large deviations from symmetry near 
b=6. It is likely that this assymetry is associated with the 
presence of a significant Galactic warp in this field. Note 
that possible evidences for an assymetry in the star counts 
has already been presented by Carney & Seitzer (1993). For 
the 2 other fields the pattern of assymetry is less clear. For 
1=180 we observe some systematic deviation near the 2 a 
level, in the range b=10 to b=30. For the field at 1=66 there 
are a few deviating points near b=4. 



3 ANALYZING THE DATA. 

3.1 Deviations from the exponential profile. 

To estimate numerically the star counts, we need to inte- 
grate the luminosity function of the stellar population and 
the density function of the disk along the line of sight. The 
luminosity function is assumed to be the same everywhere 
and is identical to the luminosity function in the solar neigh- 
borhood tabulated by Wainscoat et al. (1992). To represent 
the density profile of the disk, it is natural to look first at 
the simplest model, which is a single double exponential: 

p{R, Z) = exp [ - - M ) with R = J{x^ + y^) 

\ Rh Zh J 

It is important to note than despite its simplicity, the expo- 
nential profile seems to show en excellent consistency with 
the Z-profile of edge-on galaxies within a few scale length 
form the center (de Grijs et al. 1996). It is interesting to 
check the consistency of the data with such simple model, 
and in particular to study the shape of the deviations from 
this model which may indicate that other components are 
present. It is possible to study this model in a the full range 
of latitude covered by our data ( — 50 < b < 50, 169 line 
of sight). By probing different line of sight, with increasing 
latitudes, we probe areas of the disk more and more dis- 
tant form the Galactic plane. Thus basically, by doing this, 
we explore the Z profile of the Galactic disk. If the disk is 
really a double exponential, this Z profile should be a con- 
stant exponential, whatever the distance R from the center. 
Consequently along all the line of sight the density should 
be consistent with an exponential in Z with constant scale 
height, Hz - In case the Z profile is not exponential we should 
observe a variation of the scale height as a function of the 
latitude. To first order in Z we can develop the variation to 
the Z profile for each latitude b is equivalent to a variation 
of the scale height Hz: 

it is easy to reconstruct the function Hz{b) for the 169 line 
of sight. Along each line of sight the data have been decom- 
posed in 20 bins of magnitudes, allowing direct maximum 
likelihood determination of Hz{b) by maximum likelihood. 



It is interesting to note that the estimation of the factor 
Hz{b) along each line site can also compensate for an im- 
perfect estimation of the scale length Rh- This is obvious for 
the line of sight near I = 180, since the distance d is propor- 
tional to Z in this direction. For the line of sight near I — 240 
this is still a good approximation. An even for I = 66, one 
can check easily numerically that the factor Hz{b) can com- 
pensate with good accuracy for small discrepancies in the 
estimation of Rh- The function Hzib) has been estimated 
for the 3 fields, allowing a full reconstruction of the density 
profile of the disk. In order to make an easy comparison with 
an exponential profile, the Z profile of this density has been 
computed for different values of the distance to the center 
d. The result for the 3 fields is presented in Fig. 4,5,6. The 
most obvious difference between the 3 fields is that the ver- 
tical profile does not deviate much from the straight line 
(pure exponential) all along the line of sight in the field at 
I = 66. The situation is very different in the 2 other fields, 
where the profile becomes broader and broader with distance 
with an increasingly large deviations from an exponential. 
One could always claim that this is a bias due to the fact 
that the Z profile is constant but not exponential. However, 
we see immediately that this is not possible, since in this 
case the same discrepancy should be observed in the field 
at / = 66. There are small deviations near the top of the 
profile at / = 66, which could be due to the fact that we 
reach a distance of about 10 Kpc, for which we already ob- 
serve some flattening in the other fields. It might be also 
that in this field, source confusion affects the star counts 
in a small in a small range of latitude ( — 1.5 < 6 < 1.5). 
For the 2 other fields, the stellar density is about 10 times 
smaller, thus we do not expect any significant effects due to 
the source confusion. 

3.2 Disk models with more parameters. 

It is important to check that this thickening of the disk with 
distance that we observed is not related to some bias, due 
for instance that we did not included a thick disk. It is easy 
to check this issue by introducing a thick disk component 
in the model, and to reproduce the fitting procedure and 
the set of vertical profile in this case. Some recent determi- 
nation of the thick disk parameters indicates that the thick 
disk accounts for about 0.06 % of the local density (Buser 
et al. 1998). The scale height of the thick disk is close to 
1 Kpc while its radial scale length, is about 3 Kpc (Buser 
et al. 1998). This thick disk component was included in the 
model for the field towards the anti-center. The direction 
of the anti-center for chosen for reason of simplicity. In this 
field the amplitude of the warp is negligeable, and if the disk 
is elliptical, the effect of ellipticity is constant along the line 
of sight. The resulting set of density profiles is presented in 
Fig. 7. The addition of the thick disk component does not 
modify much the shape of the profiles. The only effect is 
only a reduction of the density in the wings of the profiles, 
but it does not affect the fiattening. Other model of thick 
disk with larger density or different scale length do not mod- 
ify the result. This can be easily understood, since due to 
its "thickness" , the thick disk becomes really important be- 
yond Z = 1.5, and most of the flattening effect we observe 
occur at Z < 1.5. Another possibility that has been used to 
reproduce the the star counts towards the anti-center is to 
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Figure 4. The vertical profile of the disk for the field at 
I = 240. The profile has been reconstructed for distances rang- 
ing from 1 to 10 Kpc along the line of sight at fe = (ac- 
cording to the limiting magnitude and the luminosity function, 
the sensitivity of the data beyond 10 Kpc is very limited). 
The corresponding distance from the Galactic center are: R = 
8.5,9.2,9.8,10.6,11.4,12.2,13.0,13,9,14.7 Kpc. The deviation from 
the exponential profile become very large beyond 13 Kpc. The 
solid lines are computed by smoothing of the data points. A line 
passing by the maxima of the different profiles has also been rep- 
resented. Note that this line is not vertical, which indicates that 
the distribution is not symmetrical, and most likely warped. 
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Figure 5. The vertical profile of the disk for the field at I = 180. 
The profile has been reconstructed for distances ranging from 1 
to 9 Kpc along the line of sight at fe = 0. The corresponding dis- 
tance from the Galactic center are: R = 9,10,11,12,13,14,15,16,17 
Kpc. The flattening of the profile with distance is obvious. The 
deviation from the exponential profile become very large beyond 
13 Kpc. A line passing by the maxima of the different profiles has 
also been represented. Note that this line is very close to verti- 
cal, which indicates shows that on the contrary to the field near 
I = 240 the warp may not be present in this field. 
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Figure 6. The vertical profile of the disk for the field at 
i = 66. The profile has been reconstructed for distances rang- 
ing from 1 to 9 Kpc along the line of sight at = 0. 
The corresponding distance from the Galactic center are: R = 
7.7,7.4,7.3,7,35,7.5,7.8,8.2,8.7,9.3 Kpc. It is important to note 
that the deviation from the exponential profile are much smaller 
in this field than in the 2 previous fields. It can be explained eas- 
ily, since in this field , even by going 10 Kpc along the line of 
sight, we do not go beyond 10 Kpc from the Galactic Center. 

used a cut-ofT of the disk to some distance of the center. For 
instance Robin et al. (1992) found that a cut-off of the disk 
at 14 Kpc was necessary in order to account for the stellar 
density in a field close to anti-center direction, near b=2.5. 
The density reconstructed with this cut-off in our model of 
the disk is presented if Fig. 8. This cut-off does not help to 
reduce the variation of profile with distance. This fact was 
easily predictable, since if we try estimate the cut-off for dif- 
ferent line of sight, the value of the cut-off is not constant 
(see Fig. 9 and 10). ft shows that a cut-off of the disk is 
not a good solution. This we can conclude that the effect of 
flattening of the profile with distance, is not a bias due to 
the simplicity of our approach, but a real effect. 

3.3 Evolution of the scale height as a function of 
the distance to the center. 

To estimate the scale height corresponding to the different 
vertical profiles an exponential has been fitted to each pro- 
file. The result for the 3 fields is presented in Fig. 11 as a 
function of the distance to the Galactic center, ft is clear 
that the 3 fields show a good consistency in the thickening 
of the profile. The effect seems to be somewhat larger at 
1=240, and is the smaller at 1=66. It is very hard to know 
what happens when going closer to the Galactic center, since 
our data do not go closer than 7.5 Kpc from the center. Un- 
fortunately the lines of sight which would allow to probe the 
disk closer to the center are not in the areas released by the 
2MASS team to date. 



3.4 Full maximum likelihood fit. 

The former analysis conducted along each line of sight has 
been particularly useful to present the data, and demon- 
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Figure 7. The vertical profile of the disk reconstructed using a 
model including a thick disk component (the field is at i = 180). 
The profile has been reconstructed for distances ranging from 1 to 
9 Kpc along the line of sight at 6 = 0. The corresponding distance 
from the Galactic center are: R = 9,10,11,12,13,14,15,16,17 Kpc. 
Note that the addition of the thick disk component has some 
influence in the range Z =1.5,2.5. However it does not change the 
shape of the profiles for Z < 1.5, and the basic features associated 
with the fiattening of the profile with distance are unchanged. 
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Figure 9. The cut-off near the Galactic along different line in 
the field at 1=180. The cut-off was estimated using a maximum 
likelihood method for different values of the disk scale length: 
3.5 Kpc (squares), 3.0 Kpc (triangles), 2.5 Kpc (crosses), and 2.0 
Kpc (stars). Note that whatever the scale length the cut-off is 
extremely variable. Thus one can conclude that a disk with cut- 
off is not appropriate to represent the stellar density along the 
line of sight. 
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Figure 8. The vertical profile of the disk reconstructed in the 
case of a model which includes a cut off of the disk density a 
R = lA Kpc. (the field is at Z = 180). The profile has been re- 
constructed for distances ranging from 1 to 9 Kpc along the line 
of sight at 6 = 0. The corresponding distance from the Galactic 
center are: R = 9,10,11,12,13,14 Kpc. Note that the effect of the 
cut-off does not help to reduce the variations of the vertical profile 
with the distance from the center. 



73 
O 



20 



o 
C 

S 16 

CO 



o 14 - 



12 



* X 



* 

X ** 

X > 
1 X X 

I A XX X 



-5 5 

latitude 

Figure 10. Same as Fig. 9 except a thick disk component (with- 
out cut-off) was included in the model. 



strate the existence of a basic feature: the thickening and 
flattening of the Z profile with distance. This effect has been 
well established, however it is also important to derive sim- 
ple analytical full maximum likelihood solution of the whole 
data set. For this full maximum solution all the line of sight 
in the range —50 < b < 50 will be fitted at the same time. 
Note that has we did previously, we will fit by maximum like- 
lihood, the decomposition in bins of magnitude (20 bins) of 
the distribution along the different line of sights, and not 



only the total counts. The effect of disk thickening will be 
investigated in 2 the fields towards the outer disk, at 1=180 
and 1=240. First, it is interesting to illustrate the fitting of 
different type of models to the data. To illustrate this issue 
we will present the result of the fitting of 4 different mod- 
els for the field at 1=180. For all the models the luminosity 
function of Wainscoat el al. (1992) has been used. The scale 
height for each population has also been taken from Wain- 
scoat el al. 1992. 
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Figure 11. The variations of the scale height of the Z profile 
as a function of the distance to the Galactic center. The crosses 
represent the field at 1=66, the stars are for the field at 1=240, 
and the triangle for the field at 1=180. 



1- ) Simple double exponential disk (Fig. 12): 
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4-) Isothermal disk + thick disk (Fig. 14): 
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5-) Disk with variable scale height (Fig. 15, 
1=180, Fig. 18, 1=240): 

/r, ^ n f R-Ro \Z-Zo{R)\\ , ,_i 

p{R, Z, I) = exp — ^ - ' ^ ,p\ " ] xil + Ar) 



Rl 



Zh{R) 



Zh{R) = 



(1 + (7? - Ro) At) X Hz R > Ro 
1 R < i?o 



6-) Disk with variable scale height and flattening 
(Fig. 17 1=180, Fig. 19, 1=240): Let's define: 



F{x,P) = 1 -exp(-/3 exp(-j;)) 



Then: 



p{R, Z, I) = exp 



(_R^^Ro\ 
\ Rh J 



X F 



Zo{R)\ 



Zht 
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(1 + At) 



F{x,P)d3 



With the same definition of Zh{R) as for the previous 
model, and: 

f l + {R~Ro) Af R> Ro 

m) = { 

I 1 R < i?o 

The best maximum likelihood fit for these 6 models are pre- 
sented in Fig. 12,13,14,15,16,17. One can note immediately 
that disk models with constant scale height do not give any 
acceptable fit of the data (Fig. 12 and 13). Even by includ- 
ing a thick disk component, it is not possible to reproduce 
the star counts (Fig. 13). By exploring different scale height 
scale length and ratio in the solar neighborhood, for the thick 
disk one cannot reproduce the star counts properly. The sit- 
uation is not improved, if instead of an exponential profile, 
one adopt an isothermal profile for the Z distribution (see 
Fig. 14, 15) However, as one can see in Fig. 16, just by adding 
a linear variation of the scale height with distance from the 
Sun, to the single exponential disk model, one can reproduce 
the star counts with much better accuracy. Note that this 
disk model with variable scale height show this good agree- 
ment with the data, with even less free parameters than the 
model with a thick disk component, which despite more pa- 
rameters is not able to reproduce the star counts. An even 
better agreement is possible if one introduces some fiatten- 
ing of the disk profile by using a "saturation" parameter /3 
(see model (6) for the definition of /3, and Fig. 16 for the fit). 
Finally, model (5) and (6) are fitted to the data in the field 
at 1=240. As one can check in Fig. 18 and 19, these simple 
model reproduce the data with quite good accuracy. Note 
that the thickening of the disk calculated using model (5) is 
the same for the 2 fields: = (1 + 0.32 {R - Ro)) x Hz- 
It is also interesting to see that the fiattening of the profile 
is larger at 1=240, something that is also visible in our pre- 
vious analysis. However something more troublesome is the 
fact that the scale length of the disk differs in the 2 fields. 
Some bias can be introduced by the presence of a spiral arm 
towards the anti-center, however the contrast in density of 
a spiral arm in the old population is small and unlikely to 
produce such discrepancies. This difference in scale heights 
is more likely due to the fact that even if our model repro- 
duces well the star counts, it is too simple to account for 
the complexity of the situation. The thickening of the disk 
is a flaring process and is likely to be non-linear, and the 
shape of the Z proflles are probably more complex than our 
simple analytical forms. To conclude with the full maximum 
fikelihood fit, one can say that the results are in very good 
agreement with the simple analysis along each line of sight 
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Figure 12. Fitting the data near 1=180 using a single disk (dou- 
ble exponential) component. The parameters of the best fit are: 
Rh = 2.45, Aw = -0.007. 



Figure 13. Fitting the data near 1=180 using a disk + a thick 
disk component. The parameters of the best fit are: Ru = 1.8, 
Aw = -0.007, Rht = 3.0, Zht = 1-0, Rs^i = 0.15 



that was conducted before. The basic features, the thicken- 
ing and flattening of the outer disk region is confirmed. 



4 DISCUSSION. 

Our first analysis, which used a first order re-construction of 
the Z profile along each the line of sight has been very useful 
to present the data and describe the basic effect. Once the 
thickening of the disk with distance had been established, 
a full maximum likelihood of the data has been performed. 
The full maximum likelihood analysis confirms the result 
obtained previously, and show that simple model, with lin- 
ear variations of the scale height and flattening are in good 
agreement with the star counts. The maximum likelihood 
analysis conflrms that the warping of the disk is present in 
the field at 1=240. As for the thickness of the disk, a lin- 
ear variation starting form the position of the Sun has been 
used to describe the warp. The maximum likelihood analysis 
shows also that a significant amount of warping is present 
in the field near 1=180, however the amplitude of the effect 
is much smaller than in the field at 1=240 (about 6 times 
smaller) . In this analysis, a large number of models has been 
tested, models with thick disk, isothermal disk models, disks 
models with cut-off..., none of them is in agreement with 
the data set, except the models including thickening of the 
disk with increasing distance from the Sun. We arrive at 
the obvious that the thickening of the disk is the only pos- 
sible explanation. Note that we did not took into account 
the metallicity effects, but since metallicity is expected to 
decrease towards the outer disk, and that stars becomes 
brighter with decreasing metallicity, the effect goes into the 
other direction. Since we observe a drop of the star counts 
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Figure 14. Single isothermal component (1=180). The parame- 
ters of the best fit are: Rh = 2.25, Aw = -0.007 



in line of sights near the Galactic plane, it cannot be caused 
by a decreasing metallicity. This possibility of a fiared and 
warped stellar disk has already been evaluated by Ewans 
el al. (1998) as a possible source of microlenses for sources 
in the LMC. This analysis demonstrates unambiguously for 
the first time the fiaring and the warping of the stellar disk. 
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Figure 15. Fitting tlie data near 1=180 using an isotliermal + 
thick disk component. The parameters of the best fit are: Rjj 
1.75, Aw = -0.007, Rht = 3.2, Zht = 1.0, Rsoi = 0.13 



Figure 16. Fitting the data near 1=180 using a disk with linearly 
variable scale height. The parameters of the best fit are: Ru = 
3.05, Aw = -0.007, At = 0.32. 



and gives an accurate description of the efltect. It is inter- 
esting to note that for other galaxies (edge-on spiral) some 
evidence of thickening of the disk vertical profile in the outer 
regions is also observed by de Grijs et al. (1996) (see Fig. 2 
in de Grijs 1996). As a whole we see that a self-consistent 
picture of the stellar disk emerge from this work: the stellar 
component is warped and flared, just like the HI component. 
The direction of the Z offset due to the warp is the same as 
in the HI maps presented by Burton & te Linkel Hekkert 
(1986). The amplitude of the effect is consistent with the HI 
to an accuracy of about 30 %. One may not forget that our 
model of the warp is very simple, in the range of distance 
investigated. Note that warp effect was treated in the linear 
approximation; it is possible that as is observed for the gas, 
the warp effect is not linear, however this linear approxi- 
mation is sufficient to reproduce the data. Concerning the 
thickening of the disk it is interesting to make a comparison 
to the fiaring of the HI disk. This comparison is presented 
in Fig. 20, scale height of the HI is taken from Wouterloot 
et al. (1990). It is clear that the amplitude of the thickening 
of the disk follow is comparable to what is observed in HI. 
This thickening of the disk, offers a natural explanation to 
the so called set of different "disk cut-off" observed for the 
Milky way (Robin et al. 1992, Ruphy et al. 1996, Freudein- 
rech 1998). due to the thickening, stars close to the plane are 
raised above the plane. This thickening gives the impression 
of a drop in the star counts at low latitudes, and fitting a 
cut-off in such density will give inconsistent results for dif- 
ferent latitudes (as illustrated in Fig. 10). This effect may 
also confuse the separation between the thick disk and disk 
component. The separation between this 2 components may 
prove very difficult. At least this work shows that clear ev- 
idences of the thick disk should be searched rather towards 
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Figure 17. Fitting the data near 1=180 using a disk with linearly 
variable scale height and flattening of the profile with distance. 
The parameters of the best fit are: Rjj = 3.05, Aw = —0.007, 
At = 0.15, Af = 0.6. 
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Figure 18. Fitting the data near 1=240 using a disk with linearly 
variable scale height. The parameters of the best fit are: Rjj = 
2.45, Aw = -0.041, At = 0.32, Zht = 1-0. 



Figure 19. Fitting the data near 1=240 using a disk with linearly 
variable scale height. The parameters of the best fit are: Ru = 
2.45, Aw = -0.041, At = 0.14, Ap = 0.81, Zht = 1-0 



the inner regions of the Galaxy. Obviously the modelisation 
of the thickening and of the warping is somewhat degener- 
ated, and it is clear that non linear thickening laws would 
fit the data also. However, in such case it is better to keep 
the simplest model that can fit the data. More accurate de- 
termination of the outer disk would require to overcome the 
degeneracy due to the unknown distances of the sources. 
Such result could be easily achieved by using contact bina- 
ries. The contact binaries are good distance indicators, and 
are found in large numbers (about 1% of the main sequence 
stars are contact binaries). These stars are ideal tools to 
probe the structure of the outer disk, and some experiment 
should be conducted in the near feature. 
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Figure 20. The distribution of the scale height of the stellar disk, 
- from first order reconstruction along the line of sight (stars) and 
from full maximum likelihood using model (5) (straight line). The 
full maximum likelihood solution as larger amplitude, this is due 
to the fact that model (5) does not include model flattening. 
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